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AIRGLOW  IMAGING  STUDIES  OF  NATURAL  AND 
ARTIFICIALLY  INDUCED  IONOSPHERIC  DISTURBANCES 

CHAPTER  1 
INTRODUCTION 

During  the  past  four  years,  a  series  of  investigations  of  "ionospheric 
variability"  have  focused  on  the  variety  of  components  that  contribute  to 
the  significant  (  20-25  )  day-to-day  variations  of  the  F-region.  These 
components  included  oeriodic  sources  (lunar  tidal  effects),  quasi -periodic 
sources  (solar  sector  boundary  effects),  and  the  essentially  random  geo¬ 
magnetic  activity  effects  (Mendillo  and  Lynch,  1979;  Schatten  and  Menriillo,  1980). 
All  of  the  above  investigations  were  carried  out  using  satellite  radio  beacon 
observations  of  total  electron  content  (TEC)  obtained  from  the  AFGL/AWS  net¬ 
work  of  TEC  monitoring  sites  (Klobuchar,  private  communication).  During  the 
final  phase  of  this  work,  our  emphasis  was  switched  to  a  different  technique 
(wide-angle  airglow  imaging)  in  order  to  investigate  two  additional  types  of 
variations  found  in  the  F-region:  (1)  the  severe  plasma  depletions  (or  "bubbles") 
that  occur  naturally  during  the  post-sunset  hours  at  equatorial  locations,  and 
(2)  the  artificially-induced  F-region  depletions  (or  "ionospheric  holes")  that 
accompany  rocket  exhaust  injections  above  200  km.  While  both  of  these  phenomena 
have  been  studied  using  radio  techniques,  airglow  imaging  diagnostics  offer  a 
new  and  complementary  approach  to  these  problems.  In  addition  to  the  usual 
benefits  of  detailed  morphology  information  associated  with  any  all-sky  photo¬ 
graphic  system,  image  processing  techniques  can  be  applied  to  airglow  images 
to  yield  additional,  quantitative  results  related  directly  to  physical  and 


chemical  processes  operating  in  the  upper  atmosphere. 

In  Chapter  2  of  this  report,  we  describe  the  low-light  level  airglow 
monitoring  system  developed  for  these  stidies.  Chapter  3  contains  a  summary 
of  the  data  reduction  and  image  processing  techniques  applied  to  the  obser¬ 
vations.  Two  experimental  campaigns  were  carried  out  to  test  the  optical 
system  for  equatorial  airglow  depletion  studies  and  rocket-induced  airglow 
enhancement  observations;  results  from  these  investigations  are  presented  in 
Chapters  4  and  5. 
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1  -K  _ 

AN  ALL  -  SKY  CAMERA  SYSTEM 

1.  Performance  Criteria 

An  all-sky  camera  system  (ASCS)  was  originally  designed  to  address  the 
needs  of  a  specific  observational  campaign:  to  monitor  airglow  depletions 
observed  in  the  equatorial  ionosphere  over  Ascension  Island.  The  airflow  data 
taken  during  this  campaign  was  to  be  used  in  conjunction  with  other  diagnostic 
techniques  to  assess  the  impact  of  equatorial  plasma  depletions  on  the  NAYSTAR 
Satellite  Navigation  System. 

The  specific  design  goals  of  the  camera  system  for  this  campaign  were: 

e 

(a)  to  record  all-sky  images  of  6300A  airglow  emissions  at  brightness  levels 

50  Rayleighs;  (b)  to  be  able  to  distinguish  features  in  the  airglow  as 

small  as  1°;  and  (c)  to  monitor  the  airglow  depletions  in  a  real  time  mode. 

Additional  design  goals  of  the  system  not  nessarily  related  to  the 

Ascension  Island  campaign  were  the  ability:  (d)  to  operate  under  adverse 

lighting  conditions  (twilight,  moon);  (e)  to  select  different  wavelength  bands 
o  o 

(e.g.,  5577A  or  7774A) ;  (f)  to  exchange  the  fisheve  lens  for  lenses  with  smaller 
fields  of  view  and  therefore  increased  resolition  1  min  of  arc);  and  (c)  to 
take  data  suitable  for  subsequent  photometric  reduction. 

2.  Design  Constraints 

The  two  most  important  constraints  in  designing  the  all-skv  camera  system 
were  cost  and  time  -  only  25  man-weeks  at  20%  effort  were  available  to  design, 
build  and  test  the  system  before  the  system  was  to  be  used  in  the  field. 


To  keep  costs  down  the  system  had  to  be  designed  using  easily  ob¬ 
tainable  "off  tlie  shelf"  components  wherever  possible.  The  decision  was 
made  to  record  the  data  photographically  rather  than  electronically  (i.e., 
television  system).  In  addition  to  reducing  the  cost,  a  photographic  approach 
has  considerably  less  bulk  and  complexity  than  a  TV  system  and  therefore  would 
be  much  easier  to  set  up  in  the  field. 

3.  Design  Details 

3  .  1  Fish-eye  lens 

Two  commonly  available  "fish-eye"  lenses  were  considered  for  the  first 
element  in  the  system.  The  N’IKKOR  8mm  FL  F/2.8  lens  produces  a  2  3mm  diameter 
image  covering  a  full  hemisphere.  The  ROKKOR-X  16mm  FI.  F/2.8  lens  for  Minolta 
cameras  has  a  40mm  diameter  image  (the  diagonal  of  35mm  frame)  also  having  a 
field  of  view  of  180°.  The  actual  choice  of  lens  was  determined  by  other 
design  constraints  in  the  system. 

3 . 2  Into r  ference  Filters 

A  major  design  problem  was  the  use  of  narrow  band  interference  filters 
with  wide  field  angle  objective  lenses.  The  use  of  narrow  band  pass  filters 
was  necessary  because  of  the  comparatively  weak  airglow  signal  ('v  100  R) 
anticipated  at  Ascension  Island,  and  a  background  continuum  of  v  3R/A  at 
6300A  (due  to  nonspectral  sources,  stars,  etc.).  A  system  would  have  to  iso¬ 
late  a  30A  wide  portion  of  the  spectrum  to  reduce  the  background  to  MOO  R. 
Filters  with  even  narrower  band  passes  than  this  would  be  necessary  if  the 
system  had  to  perform  under  adverse  lighting  conditions  (e.g.,  scattered 
moonlight  could  produce  't  20R/A) . 


5 


The  difficulty  in  using  interference  filters  to  limit  the  band  pass  of 
an  all-sky  system  is  that  the  band  pass  of  an  interference  filter  shifts  to 
shorter  wavelengths  for  rays  that  are  skew  to  the  filter's  normal  (see  figure 
1) .  The  severity  of  this  shift  depends  on  the  index  of  refraction  of  the 
materials  used  to  construct  the  filter  (see  figure  2). 

Two  possibilities  for  using  interference  filters  with  all-sky  imaging 
systems  were  investigated.  Figures  3  and  4  show  the  two  designs  considered 
in  schematic  form. 

In  Figure  3,  the  filter  is  placed  in  the  diverging  rays  from  the  fish- 

eye  lens  after  the  >-ays  have  passed  through  the  focus.  A  field  lens  has  been 

added  at  the  image  plane  to  make  the  central  axis  of  the  diverging  cones  of 

light  parallel  to  the  optical  axis  of  the  system.  This  approach  is  similar 

to  that  of  Mende  et  al.  (1976).  The  maximum  angle  that  any  ray  makes  with 

the  filter  normal  in  this  svstem  is:  ex  =  TAN  1 ( — )  .  For  either  of  the  lenses 

2F 

considered  above  a  -  10  .  The  sizes  of  the  filters  using  this  approach  would 

be  approximately  the  same  as  the  sizes  of  the  images  formed  by  the  fisheve 

lenses.  Referring  to  Figure  5,  it  can  be  seen  that  the  filter  passband  would 

o 

have  to  be  <v>  40  A  and  centered  at  6310A,  so  that  all  the  rays  will  pass  within 
the  907  (normalized)  transmission  region  of  the  filter.  Note  that  in  this 
design  the  rays  from  the  edge  of  the  field  of  view  of  the  fisheye  are  treated 
the  same  as  those  from  the  center  of  the  field  of  view  (see  Fig. 3). 

Figure  4  shows  an  alternate  method  of  treating  the  problem.  In  this 
approach,  the  diverging  rays  from  the  image  are  collimated  before  passing  through 
the  filter.  To  ensure  that  all  the  light  from  the  fisheye  enters  the  collimator 
lens,  a  field  lens  is  again  employed.  The  focal  length  of  the  field  lens  is 
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such  that  it  images  the  aperture  stop  of  the  fisheye  onto  the  filter.  The 
filter  sees  a  bundle  of  rays  from  the  center  of  the  field  of  view  all  parallel 
to  the  optical  axis.  For  rays  originating  near  the  edge  of  the  field  of  view, 
the  filter  sees  a  parallel  bundle  of  rays  whose  axis  makes  an  angle  to  the 
optical  axis  equal  to  one-half  the  angular  size  of  the  image  as  viewed  from 
the  collimator.  It  is  obvious  that  this  angle  can  be  made  arbitrarily  small 
by  placing  the  collimator  far  away  from  the  image.  This  angle  is  also  mini¬ 
mized  by  starting  with  a  fisheye  lens  that  produces  a  small  image.  In  prac¬ 
tice,  however,  the  collimator  can  not  be  placed  at  great  distances  from  the 
image  because  the  diverging  beam  from  the  image  gets  too  large.  Therefore, 
the  practical  limit  on  the  focal  length  of  the  collimator  is  set  by  how 
large  a  filter  can  be  obtained.  Without  going  through  the  extraordinary 
procedure  of  making  a  mosaic  filter,  the  practical  upper  limit  for  interfer¬ 
ence  filters  with  't  6A  passbands  was  found  to  be  't  100mm  in  diameter  (the  main 
constraints  being  cost  and  quality  control).  The  selection  of  this  filter  size 
and  the  NIKKOR  Fisheye  lens  defined  the  minimum  passband  possible  for  a  filter. 
The  maximum  spacing  between  the  image  formed  by  the  fisheye  lens  and  a  100mm 
diameter  collimator  is  285  mm.  To  separate  them  more  than  this  would  allow  the 
diverging  beam  to  overfill  the  collimator.  With  a  collimator  focal  length  of 
285  mm  and  an  image  radius  of  11.5  mm,  the  maximum  angle  anv  ray  would  make 
with  the  filter  normal  is  2.3  .  Referring  again  to  Figure  2,  this  angle  of 
incidence  would  result  in  a  wavelength  shift  of  approximately  1.5A.  To  ensure 
that  all  of  the  rays  pass  through  the  filter  within  the  907  (normalized) 

transmission  region,  the  FWHP  passband  of  the  filter  would  have  to  he  at  least 

o 

A301A  (sec  figure  6).  The  center,  wave  1  ength  of  these 


3A  and  be  centered  at 
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filters  shift  to  longer  wavelengths  as  temperature  is  increased  with  a 
coefficient  of  'v  0.2A/°C.  Therefore,  the  filters  are  housed  in  a  temperature 
controlled  oven. 

The  closest  match  to  the  ideal  collimator  lens  for  the  system  was  an 

ACHROMAT  88mm  in  diameter,  24Lmm  focal  length.  This  slightly  shorter  focal 

o 

length  increased  the  maximum  angle  of  incidence  at  the  filter  to  2.7  and 
therefore,  the  minimum  pass band  of  the  filters  to  be  used  with  the  system 
increased  to  6A. 

At  ter  passing  through  the  filter,  the  collimated  light  is  rc-i.ma  «ed  V  v 
.i  Jena  identical  to  the  e  llimator  described  above  .  This  filtered  image  !  s 
tin-  same  :;:::e  a.;  the  imam  termed  by  tin  fi shove  lens  and  is  slightly  concave 
toward  the  «>'■•  jest  ive  1  en  -  .  pro  Mom  oas  j  1  \  !■  sailed  by  using  an  appropriate!  v 

curved  fih>r  ,,t  i  race  plate  set  in  c'!  against  the  (lane  surface  of  the  input 

of  t  hi'  image  intorsi  :  i  vr  .  r  ibed  :  n  t  :  •.  f-i  '  low  in;.',  section. 

3.3  imre  e  int.  ns.  i:  ier 

The  r  i  ‘un  i  raiioii  t  that  (he  s\  steu.  fo-.o'.-a  !  ‘  and  the  fact  t  ••at  the  air glow 

feat  ire  •  to  he  -t'.idi'-d  >ro  known  :  ■  move  at  I'/min  at  tin  :;cn  ith  (Weber ,  et  a  t 

!"?>())  limited  the  expo  r-  time  >  o  *  ••  -  :  ban  up  mir.uto . 

This  c>ir  si  ra  ir.t ,  coupled  with  t  he  expected  bright  ness  of  the  equatorial 
uirv low  (jl)  -  l DOR) ,  mandated  the  use  of  an  image  intensifying  device  to  get 
usable  t  ilm  densities.  two  basic  r  pes  of  image  intonsif iers  were  considered 
for  the  camera  system:  A  23mm  3-stage  so-called  "first  generation"  image  tube 
and  a  second  generation  m  icrochanne  l  plate  (MCP)  intensif ier .  Both  tubes 
have  similar  gain,  noise,  and  resolution  specifications.  However,  the  23mm 
MCP  intensifier  (military  specification:  MX9644/UV)  is  one-third  the  length 
of  the  three  stage  tube,  has  an  integral  H.V.  power  supply  (needing  only 


external  2.5V  battery  for  power)  and,  most  importantly,  has  antiblooming 
characteristics  that  make  it  particularly  well-suited  for  use  in  a  situation 
where  a  large  dynamic  range  might  be  encountered (a . g . ,  the  moon's  image 
simultaneously  with  dim  airglow) .  Figure  7  shows,  in  schematic  form,  the 
operation  of  a  second  generation  image  intensifier. 

Photons  enter  the  fiber  optic  faceplate  and  strike  the  S20/25  photocathode 
which  is  deposited  on  the  inside  surface  of  the  faceplate.  Photoelectrons 
liberated  at  the  cathode  are  then  accelerated  and  electrostatically  focussed 
on  the  microchannel  plate  (MCP)  electron  amplifier.  The  MCP  consists  of  a 
parallel  array  of  hollow  glass  cylinders  about  10pm  in  diameter  and  about  1  mm, 
or  less,  in  length.  The  inside  walls  of  the  cylinders  are  coated  with  a  secondary 
emitting  material.  When  a  primary  electron  enters  a  cylinder,  it  cascades  down 
the  channel  producing  additional  secondary  electrons. 

Bonded  to  the  output  end  of  the  MCP  is  a  thin  aluminum  barrier  and  a 
P1/P39  phosphor  screen  which  converts  the  electron  current  back  into  photons. 

Data  sheets  describing  the  characteristics  of  the  intensifier  used  in  the 
all-sky  system  are  included  as  an  Appendix  to  this  Chapter. 

3 .A  Camera 

To  photographically  record  the  image  from  the  output  phosphor  of  the 
intensifier,  a  transfer  lens  system  and  a  standard  35mm  single  lens  reflex 
camera  were  used.  A  more  efficient  method  of  coupling  this  image  to  the  film 
would  be  to  have  the  emulsion  of  the  film  pressed  against  the  fiber  optic 
output  of  the  image  tube.  This  method  permits  all  of  the  light  from  the  phosphor 
to  enter  the  emulsion.  However,  a  film  transport  mechanism  to  accomplish  this 
is  a  custom  made  device  -  and  most  importantly,  if  such  a  system  is  employed, 
the  image  cannot  be  viewed  in  real  time. 
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The  penalty  for  optically  transferring  the  output  phosphor  to  the  film 
is  rather  severe  -  only  about  10%  of  the  light  from  the  intensifier  can  be 
coupled  to  the  film  emulsion  with  transfer  lens  systems  operating  near  F/l. 
Special  "fast"  lenses  have  been  designed  to  operate  at  conjugate  ratios  of 
1:1  to  address  this  demanding  situation.  However,  it  was  found  that  com¬ 
monly  available  "fast"  lenses  for  35mm  cameras  could  be  used  in  tandem  pairs 
to  produce  as  efficient  transfer  systems.  In  such  a  system  the  lenses  are 
placed  front  to  front  so  that  the  lenses  are  used  at  their  designed  infinite 
conjugate  ratios.  The  first  lens  acts  as  a  collimator  while  the  second 
reimages  the  collimated  beam  on  to  the  film.  Such  a  system  suffers  from 
vignetting  but  with  care  this  can  be  minimized.  In  the  current  system  an 
85mm  F/l. 8  Nikon  lens  is  used  as  the  collimator  and  a  55mm  F/l. 2  Nikon  lens 
is  used  on  the  camera.  This  combination  of  focal  lengths  produces  an  image 
on  the  film  2/3  the  size  of  the  image  on  the  output  of  the  intensifier  and 
keeps  the  total  vignetting  of  the  system  to  less  than  a  factor  of  3  over  the 
image  (see  Figure  8) .  When  used  in  the  all-sky  mode  with  the  camera  pointed 
at  the  zenith,  this  vignetting  is  partially  offset  by  the  Van  Rijn  brightening 
at  large  zenith  distance.  However,  when  the  system  is  used  with  a  narrow 
(-60°)  field  of  view  lens  the  vignetting  is  very  noticeable. 

A  Nikon  F  camera  body  with  a  motor  driven  250  exposure  back  is  used  to 
record  the  images.  The  camera  back  has  been  modified  to  allow  data  (time, 
date,  etc.)  to  be  recorded  along  one  side  of  each  exposure.  The  camera  is 
operated  by  an  external  intervalometer  built  specifically  for  the  Ascension 
Island  campaign.  The  interval  between  exposures  can  be  set  to  2’^,  5  or  10 
minutes  and  the  exposure  time  can  be  varied  by  factors  of  2  from  1/1000  sec 


FIGURE  8 
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to  32  sec.  A  schematic  of  the  completed  system  and  the  system  set  up  in  the 
field  are  shown  in  figures  9  and  10. 

4.  Measured  Performance 

4 . 1  Sensitivity 

Sensitivity  is  defined  here  as  the  minimum  exposure  in  Rayleigh-seconds 
(RSEC)  to  produce  a  density  0.3  above  the  combined  density  of  the  base  of  the 
film  and  fog  and  tube  noise. 

Several  film  were  investigated  until  KODAK  TRI-X  developed  in  D-19  for 
5  minutes  was  selected  as  best  for  overall  speed,  resolution,  and  availability. 
Using  this  film  and  the  above  definition  the  sensitivity  of  the  system  is 
800-1000  RSEC. 

4 . 2  Resolution 

Field  tests  of  the  system  have  demonstrated  that  the  system  can  resolve 
features  smaller  than  1°  at  the  center  of  the  field  of  view  of  the  f isheve  lens. 
This  corresponds  to  about  12  lp/mm  at  the  image  plane  of  the  system. 

4 . 3  Real  time  viewing  of  weak  airglow 

At  Ascension  Island  airglow  depletions  could  clearly  be  seen  in  the  view¬ 
finder  of  the  camera.  A  clear  focusing  screen  and  a  6x  magnifying  viewer  were 
used  to  maximize  the  brightness  of  the  image.  An  RCA  ultracon  T.V.  camera 
"looking"  into  the  camera  viewfinder  proved  not  to  have  the  sensitivity  to 
produce  usable  images. 

4 . 4  Operation  under  high  background  conditions 

During  the  Ascension  Island  campaign  no  adverse  lighting  conditions  were 


encountered.  The  campaign  was  scheduled  so  that  the  moon  was  below  the  horizon 
during  the  hours  of  interest.  However  the  narrow  bandwidth  of  the  filters  used 
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(12A)  allowed  useful  airglow  measurements  to  be  taken  starting  at  solar 
zenith  angle  'U  100°. 

The  system  was  used  to  monitor  a  rocket  launch  for  Vandenberg  AFB  on  June 
23,  1981.  (see  Chapter  5),  This  was  an  "experiment  of  opportunity"  and  as  such 
the  launch  data  and  time  were  dictated  by  the  needs  of  the  spacecraft  orbit  and 
not  uiir  desire  to  monitor  the  ionospheric  effects  of  the  launch  vehicle.  The 
3rd  quarter  moon  v>a s  (>0  away  from  the  airglow  cloud  and  was  imaged  onto  the 
photocat  1  '.it*  f  the  image  intensifier  by  the  180°  f i.sheye  lens.  An  •  ample 
of  a  photograph  taken  during  this  event  is  hown  in  Figure  11.  It  car.  he  see:', 
that  the  moon's  image  did  not  prevent  useiul  data  from  being  taken. 

■1.5  T nte r c hangeabl  e  _  lenses 

For  cite  June  rocket  campaign  the  system  was  configured  with  a  60°  iield 
view  objective  lens.  8  minutes  or  so  after  the  launch  the  airglow  cloud 
exceeded  this  field  and  the  *50°  fie.d  len-  wu  replaced  by  the  f  isheye  —  a 
pr.u  ■.  ss  that  takes  about  5  seconds.  A  series  of  photographs  taken  with  the 
•>c°  f  it: Id  of  view  lens  are  shown  Ln  Chapter  5. 

4 . 0  Photometric  a re  u  racy  of  the  system 

Chapter  3  discusses  at  length  the  photometric  reduction  techniques  applied 
to  the  data  taken  with  the  intensified  camera  system. 

5 .  Fu t  u r e  Mod i f i ca  t  ion s 

Some  mod  i  f  ir  it  ion-,  planned  to  t  he  intensified  svstem  include  : 

1.  addition  of  difterent  wavelength  filters; 

2.  increasing  the  r  'upline  efficiency  between  the  output  of  the  inten- 
sif  ier  and  t  lie  film; 

3.  the  addition  of  i  more  sensitive  TV  camera  to  the  view  tinder; 

4.  a  more  flexible  int erva 1  one t er  ; 

5.  addition  of  a  thermostat  h  .il  lv  controlled  over,  tor  the  filter  wheel. 
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CHAPTER  3 

DATA  REDUCTION  AND  IMAGE  PROCESSING 

1 .  Int  rodue  t ion 

In  this  chapter,  a  discussion  of  the  various  types  of  data  reduction  and 
analysis  will  be  given  for  the  data  gathered  during  the  Atlas-F  rocket  launch 
of  June  1981,  and  the  Ascension  Island  equatorial  airglow  depletion  campaign 
of  January  1981.  The  nature  of  the  data  reduction  concerned  the  processing  of 
data  captured  on  >3  rm.  film  to  same  foim  of  visual  display  capable  of  relating 
qualitative,  as  well  as  quant  itat iv--  in  format  ion.  Initially,  prints  were  made 
of  the  frames  describing  the  events  and  simple  overlays  were  created  to  charac¬ 
terize  the  gross  features.  I.ater,  these  same,  data  were  digitized,  and  through 
a  library  of  computer  programs,  were  produced  graphically  in  a  manner  more 
conducive  to  scientific  investigation. 

The  first  method  of  analysis  for  the  airglow  .images  was  developed  for  the 
Ascension  Island  data.  A  series  of  three  computer  generated  overlays  were 
developed  for  use  on  the  photographic  prints  or  the  airglow  depletions.  The 
first  overlay,  created  from  a  starfield  taken  at  a  known  time,  was  used  to 
determine  the  altitude  and  azimuth  of  any  given  point  in  the  f ie Id -of -v iew . 

The  second  overlay,  computed  for  a  specific  airglow  height,  could  be  used  to 
determine  the  subionospher ie  geographic  coordinates  of  a  depletion.  Finally, 
the  last  overlay  produced  a  set  of  magnetic  meridians,  again,  at  a  given  height 
above  the  Earth.  These  overlays  could  be  used  to  give  quantitative  information, 
such  as  airglow  depletion  drifts,  spatial  extents,  and  alignments  with  respect 
to  geomagnetic  coordinates.  In  addition,  for  the  rocket  data  gathered,  one 
similar  type  of  overlay  was  devised.  This  overlay  displayed  the  trajectory  of 
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the  Atlas-F  rocket,  projected  onto  an  altitude-azimuth  grid  with  times- 
of-flight  and  corresponding  rocket  heights. 

2.  Digital  Analysis  Methods 

Digitization  of  ’’he  data  by  the  Eikonixscan  system  yielded  3  12  x  312 
arrays  of  8  bit  transmission  numbers.  It  was  decided  to  reduce  the  size 
of  these  large  arrays  for  several  reasons.  Primarily,  the  resol  at  lor  the 
all-sky  monitoring  system  with  the  fish  eye  lens  (180  degrees;  ranged  from 
1.0  to  2.0  degrees,  and  about  three  times  greater  with  t ae  narrow  field  lens 
'  nO  degrees.'.  Thus ,  an  arras-  "12  s  3 1 2  in  size  would  ove  r  sample  each  frame, 
in  effect  creating  detail  that  -was  beyond  t  he  opt  ical  limit  ol  resoiut  i  <  •» . 

In.  addition,  Lite  redic  Lion  o'  tile  array  served  to  smooth  out  any  anomalous 
pixi  !  to  pixel  d  i  si  .repane  im  and/or  sc  r. at  cites  or  emulsion  detects  (e.z., 

.  Lnheies,  ittist  ,  etc.)  on.  t  In  i  i  Im  it  sel :  .  Finally,  from  a  practical  stand¬ 
point  a  smaller  array  was  much  faster  to  process  in  digital  computations. 

In  order  to  reduce  the  array  size,  and  still  maximize  the  three  criter¬ 
ion  mentioned  above,  i.e.,  (1)  to  approach  a  size  comparable  to  t he  resolution 
without  creating  or  destroying  detail,  (2)  to  smooth  our  troublesome  pixel 
noise  without  suppressing  meaningful  data,  and  (3)  to  choose  the  smallest 
possible  array  to  promote  speed  and  efficiency  in  processing,  it  was  decided 
to  use  a  1 28  x  128  pixel  array.  This  size  was  obtained  by  a  simple  averaging 
of  adjacent  4x4  regions  to  create  the  final  128  x  128  format  used  in  all 
following  discussions. 

In  addition  to  the  data  taken  in  the  field,  control  data  were  also  digitized 
so  that  a  means  to  go  from  transmission  units  to  intensity  units  could  he 
established.  Using  the  same  all-sky  camera  set-up,  a  standardized  light  source* 
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was  photographed  across  its  range  of  graduated  intensities  using  different 
exposure  times  in  order  to  simulate  the  conditions  under  which  the  field  data 
were  gathered.  After  digitization  of  these  uniformly  illuminated  frames, 
both  with  and  without  a  1.0  neutral  density  filter,  the  center  5x5  pixel 
region  of  each  array  was  averaged  (to  avoid  possible  vignetting  effects). 

The  end  product  of  such  a  procedure  was  a  transmission-to-intensity  look-up 
table  for  different  exposure  times  and  filter  combinations.  To  find  the 
intensity  corresponding  to  a  given  transmission  number,  it  was  necessary  to 
account  for  whether  the  frame  was  digitized  witli  or  without  a  filter,  and  what 
the  corresponding  exposure  time  was;  the  table  then  yielded  the  correct  inten¬ 
sity  value. 

Before  the  intensity  transformation  could  be  carried  out,  it  was  necessarv 
to  correct  for  any  non-uniformities  present  in  the  digitizing  system.  This  was 
accomplished  by  digitizing  a  field  with  no  negative  in  the  film  carrier.  This 
so-called  "flat-field"  array,  in  the  absence  of  non-uniformities  of  illumination 
and  diode  sensitivities,  should  have  the  same  value  everywhere.  The  transmission 
numbers  in  the  actual  "flat-field"  array  showed  about  15'  variation  over  t hi' 
field  digitized.  By  normalizing  this  flat  field  frame  to  the  highest  transmission 
value,  an  array  was  obtained,  which  when  divided  into  any  other  data  arrav  in 
transmission  units,  would  correct  for  uneveness  of  illumination  oi  the  negative 
being  d igitized . 

In  addition  to  the  "flat-field"  correction,  it  was  clear  from  the  design 
of  the  all-sky  imaging  system  that  internal  effects,  such  as  t he  vignetting 
encountered  in  reimaging  the  output  phosphor  of  the  intensit  jer,  as  well  as 
external  factors,  such  as  the  van  Rhijn  effect,  must  he  addressed.  hue  to 
experimental  circumstances,  these  problems  were  approached  in  different  wavs 
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m  regards  to  the  Atlas-K  rocket  launch  and  to  the  Ascension  Island  depletion 
campaign.  First,  tae.  rocket  launch  data  analysis  for  vignetting  correction 
is  discussed. 

l'o  produce  a  vignetting  array,  a  frame  was  digitized  that  contained  a 
uniform  signal,  i.e.,  a  frame  which  was  exposed  while  the  camera  system  was 
monitoring  an  evenly  illuminated  diffuse  source.  While  the  array  was  in 
transmission  units,  the  standard  "flat-field"  correction  was  applied,  and, 
as  described  net ore,  the  appropriate  transformation  to  intensity  was  per- 
torrued.  Again,  as  with  the  "flat-field"  array,  this  the  so-called  "vignetting" 
arrav,  was  normalized  to  its  highest  value.  Similarly,  it  was  then  used  to 
correct  other  data  arrays  at  ter  they  were  transiormed  to  intensity  units. 

Be  lore  tiifs  correction  was  actually  applied,  it  became  apparent  that  the 
vignetting  array  was  a  bit  noisy,  though  it  did  resembLe  closelv  the  characteris¬ 
tic  buils-eye  pattern  ot  an  analytically  defined  vignetting  function.  After 
repeated  attempts  at  two-dimensional  spline  tits  it  was  decided  to  use  a  simple 
i  x  j  running  average  scheme.  This  type  of  averaging  involved  moving  one 
pixel  at  a  time  and  averaging  the  surrounding  lb  pixels.  This  method  was  by 
tar  the  most  successful  at  smoothing  large  spikes,  while  at  the  same  time,  it 

was  equally  ettective  in  maintaining  detail  and  a  close  resemblance  to  the  raw  ». 

uusinoothed  aaLa.  Thus,  tnis  method  resulted  in  creating  a  workable  vignetting 
array  wn  it:  h  was  an  essential  tool  tor  contouring  routines. 

For  the  Ascension  island  airglow  depletion  photographs ,  different  physical 
conditions  led  to  a  revised  analysis  of  vignetting  and  van  Khijn  effects.  In 
contrast  to  the  rocket  data,  the  Ascension  Island  data  were  taken  onlv  with  the 
IdO  degree  t ish-eye  lens.  A  new  approach  was  taken  to  correct  for  vignetting 
and  van  Khijn  effects  together. 


Inspection  of  all  the  photographs  digitized  showed  that  no  one  area  of 
the  sky  was  preferred  by  depletion  activity.  After  choosing  two  representative 
days,  (Feb.  1  and  2,  1981),  each  frame  was  corrected  for  "flat-field"  and 
put  into  intensity  units.  Then,  all  the  image  arrays  for  a  night  were 
averaged  together  and,  finally,  the  two  night's  averages  were  combined  in  a 
grand  average,  the  so-called  "average  frame."  This  average  frame  was  then 
smoothed  by  the  5x5  pixel  averaging  method  described  above.  The  end  result 
was  an  "average  frame"  that  yielded  the  combined  mean  total  of  the  vignetting 
of  the  system,  the  van  Rhijn  effect  of  the  atmosphere,  the  background  intensity 
of  the  airglow,  and  the  stationary  objects  which  were  unwanted,  such  as  tele¬ 
phone  wires  (see  figure  1). 

The  average  frame  was  used  in  two  separate  ways.  In  the  first  technique 
the  average  frame  was  subtracted  from  the  data  array  to  form  a  "difference 
image."  In  such  a  format,  all  stationary  objects  and  constant  airglow  features 
went  to  zero,  and  if  a  constant  factor  was  added  to  the  array,  the  depletions 
appeared  as  indentations  along  a  continuum.  While  this  method  did  not  really 
address  the  issues  of  vignetting  and  the  van  Rhijn  effects,  nor  was  it  capable 
of  producing  absolute  intensities,  it  was  quite  effective  when  used  to  display 
the  essential  characteristics  of  airglow  depletions  as  departures  from  mean 
conditions . 

The  second  technique  took  the  average  frame,  normalized  it  to  the  peak 
value  in  the  array,  and  divided  it  into  each  of  the  data  arrays.  This  resulted 
in  a  good  correction  for  the  combined  effects  of  vignetting,  van  Rhijn,  and 
local  time  difference  from  east  to  west,  but  did  not  totally  "erase"  the  pre¬ 
sence  of  the  stationary  objects  in  the  f ield-of -view.  This  method  was  essential 
for  determining  the  true  slopes  of  the  walls  of  the  depletions  and  was  the 


33 


preferred  correcting  function  to  apply  for  contouring  purposes. 

In  summary,  for  both  the  Atlas-F  and  the  Ascension  Island  projects, 

the  data  reduction  followed  along  quite  similar  lines  (see  figure  2)  .  Each 

photograph  was  digitized  to  a  512  x  512  pixel  array  that  reduced  subsequently 

to  a  128  x  128  array  by  a  straight  4x4  averaging  in  transmission  units. 

These  arrays  were  then  corrected  for  flat-field,  converted  to  intensities  via 

an  appropriate  look-up  table,  and  then  corrected  for  the  corresponding  vig- 

o 

netting  and  van  Rhijn  function.  By  inspection  of  control  filters  at  6200A 

for  both  events,  it  was  decided  that  any  background  light,  aside  from  the 
o 

6300A  airglow,  was  below  the  sensitivity  of  the  system  and  could  therefore 
be  ignored  in  the  reductions.  The  final  data  arrays  were  usually  still  some¬ 
what  noisy  to  be  of  practical  use  in  further  computations.  Thus,  in  nearly  all 
cases,  a  5  x  5  running  average  technique  was  applied  to  help  smooth  the  data 
for  final  analysis  or  display. 

3.  Graphical  Display  of  the  Final  Data  Arrays 
The  ionospheric  hole  created  by  the  rocket  was  the  first  data  to  be 
graphed  and  had  a  few  problems  associated  with  it.  The  airglow  burst  recorded 
was  on  the  order  of  14  kilorayleighs ,  far  above  an  almost  negligible  background. 
In  order  to  increase  the  effective  dynamic  range  of  the  scanning  system,  these 
frames  were  digitized  both  with  and  without  a  filter.  This  was  done  in  the 
hope  that  a  linear  relation  could  be  found  between  the  two  resulting  arrays. 
However,  due  probably  to  scattered  light  in  the  digitizing  system,  no  such 
relation  was  clearly  defined.  This  resulted  in  the  high  end  of  the  intensity 
scale  (i.e.,  the  densest  part  of  the  negative)  to  exhibit  a  fair  amount  of 


unce  r ta inty . 


Since  a  calibrated  photometer  was  used  in  conjunction  with  the  imaging 
data  taken  during  the  rocket  experiment,  absolute  values  of  intensities  could 
be  fixed  to  the  data.  Simply  by  knowing  where  the  photometer  was  pointing, 
the  corresponding  pixel  in  the  array  could  be  monitored  over  time  and  essentially 
matched  to  the  photometer  tracing  along  the  same  interval  to  produce  a  means 
of  calibration  to  absolute  intensity.  These  intensities  were  then  transformed 
to  log  ( ir.tensit  ies)  due  to  the  wide  range  of  values  encountered  during  the 
event . 

In  graphing  these  data,  it  became  necessary  to  suppress  extraneous  signals 
from  the  intensifier  that  were  bevond  the  edge  of  the  actual  f ield-of-v iew  of 
the  camera.  The  first  contouring  method  involved  drawing  isophotes  at  certain 
prescribed  levels  on  a  two-dimensional  graph  ('x'  pixel  vs.  'v'  pixel).  This 
tvpe  of  plot  showed  well  the  wav  in  which  the  hole  progressed  in  both  shape  and 
size  over  time.  Specifically,  for  example,  a  certain  intensity  level  could  he 
watched  to  estimate  expansion  rates.  The  other  tvpe  of  contour  graph  consisted 
of  taking  this  two-dimensional  grid,  projecting  it  in  the  "z"  direction,  and 
then  dropping  a  net  or  fine  mesh  atop  it.  Tn  this  wav,  a  three-dimensional 
representation  of  the  phenomenon  could  be  presented,  where  the  hole  appeared 
as  n  lull  above  the  background  airglow.  The  best  leature  of  this  technique  was 
that  the  graph  could  he  rotated  to  allow  various  viewing  angles  so  that  details 
in  the  structure  of  the  hole,  i.e.,  t  tie  observed  "horshoe  pattern",  could  be 
t ol lowed  as  it  changed  over  time.  Plots  of  all  the  digitized  Mi  degree  frames 
are  given  in  Figure  i. 

Many  of  the  digital  display  schemes  developed  for  the  rocket  launch  effects 
were  subsequent lv  used  for  the  display  of  t  tie  airglow  depletions  recorded  at 


Ascension  Island.  Because  the  ranee  of  intensities  of  the  depletions  were 
much  less,  no  filter  was  needed  in  digitizing  and  thus  coni  idence  in  intensity 
values  were  good.  However,  without  the  benefit  of  photometer  readings,  the 
best  that  could  be  arrived  at  were  relative  intensity  levels.  For  the  most 
part,  as  mentioned  before,  the  use  of  log  (intensities)  was  not  required  as 
the  range  of  values  was  not  as  large  as  the  rocket  event. 

As  with  the  rocket  data,  both  types  of  contouring  routines  were  used  for 
the  equatorial  depletion  data.  Again,  the  signals  exterior  to  the  actual 
; i e ld-of -view  'were  suppressed  and  viewing  angles,  where  necessary,  could  be 
altered  to  sight  down  depletion  "valleys."  In  the  end,  each  type  of  presenta¬ 
tion  was  useful  to  show  the  finer  structures  of  the  features  and  to  show  the 
temporal  progression  of  the  depletions. 

In  addition  to  these  contour  plots,  a  third  type  of  graph  was  developed 
for  the  Ascension  Island  data.  Before  performing  the  5x5  running  average 
on  the  corrected  arrays,  the  line  of  numbers  running  from  the  east  horizon  to 
the  west  horizon  through  the  zenith,  the  so-called  "prime  vertical"  was  extracted. 
I’sing  the  fact  that  Ascension  Island  is  close  to  the  equator  and  that  there 
existed  a  linear  relation  between  pixel  number  and  an  are-d  istanee  iron  the 
zenith,  it  was  possible  to  convert  pixel  number  to  distance  at  any  given  height 
above  the  Earth.  These  prime  vertical  plots  (relative  intensity  vs.  distance) 
then  could  be  used  to  calculate  drift  velocities  ns  well  as  the  slopes  el  the 
depletion  walls. 
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FIGURE  2.  Imaging  Processing  Technique 
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FIGURE  3.  Digital  image  processing  of  airglow  photographs 
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CHAPTER  4 

AIRCLOW  CHARACTERISTICS  OF  EQUATORIAL 
PLASMA  DEPLETIONS 


I.  INTRODUCTION 

Airqlow  observations  have  long  been  used  to  study  the  structure  and 

dynamics  of  the  upper  atmosphere  (Cabannes,  1935;  Bates,  1946).  The  oxygen 

o 

"red  line"  emission  at  6300  A  offers  a  particularly  convenient  way  of 
monitoring  processes  in  the  lower  portions  of  the  ionospheric  F-region 

(Chamberlain,  1961;  Peterson  et  al.,  1966).  Hays  et  al.,  (1978)  have  pro- 

c 

vided  a  recent  description  of  the  various  sources  and  sinks  of  6300  A 
radiation  in  the  ionosphere.  During  the  nighttime  hours ,  photochemical 
equilibrium  exists  between  the  production  of  0(1D)  atoms  via  the  F-region 
loss  process 

o+  +  o:  I  o/  +  0  (1) 

0:+  +  e-  “  k0( 1 D)  +  (2-k)0  (2) 

and  the  decay  of  0 ( 1 D )  by  emission  and  quenching.  The  charge  transfer  re¬ 
action  (■>•),  the  dissociative  recombination  rate  (u),  and  the  efficiency 
factor  (k)  for  the  number  of  0(!D)  states  created  per  plasma  recombi nation 

are  considered  to  be  well-known  parameters  (Hays  et  al . ,  1978),  and  thus 

o 

the  overall  emission  at  6300  A  depends  on  the  height  profiles  of  0:  and  e'. 
Since  the  equatorial  ionosphere  exhibits  a  well-known  "anomalous  pattern" 
for  the  height  (hmax)  and  magnitude  (Nmax)  of  its  peak  density,  while  the 

neutral  0-,  profiles  are  relatively  uniform  within  ±  15°  dip  latitude  of  the 

o 

magnetic  equator,  the  airglow  emission  at  6300  A  depends  almost  exclusively 
upon  the  instantaneous  structure  of  the  Ne(h)  profiles.  The  morphology  of 
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the  F-region's  equatorial  anomaly  has  been  reviewed  recently  by  Anderson 
(1981)  and  Anderson  and  Roble  (1981);  their  results  show  that  during  pre¬ 
midnight  hours,  the  crests  of  the  equatorial  anomaly  have  hm  values 
near  300  km,  while  at  the  equator  h  aDoroaches  6C0km.  This  extra- 
ordinary  height  variation  accounts  for  the  prominent  intertropical  arcs 

O 

or  bands  of  6300  A  airglow  associated  with  the  anomaly,  as  first  discussed 
in  detail  by  Babier  et  al.  (1961). 

In  recent  years,  airglow  studies  of  the  equatorial  ionosphere  have 
concentrated  on  the  departures  from  uniform  airglow  contained  within  the 
broad  intertropical  arcs.  This  new  approach  resulted,  in  part,  from  a 
new  observing  technique  developed  by  Mende  and  Father  (1976)  and  Mende 
et  al.  (1977)  in  which  a  TV  system  operating  in  an  integrating  mode,  with 
a  wide  field  of  view  and  a  narrow  spectral  bandwidth,  records  monochromati c 
all-sky  observation  of  auroral  and  airglow  emissions.  Using  this  type  of 
instrument  installed  in  the  Air  Force  Geophysics  Laboratory's  Airborne 
Ionospheric  Obeservatory ,  Weber  et  al .  (1978,  1980)  obtained  a  remarkable 
series  of  "equatorial  airglow  depletions"  from  all-sky  imaging  observations 

O  0 

at  6300  A.  The  regions  of  decreased  6300  A  intensity  were  north-south 
aligned,  had  east-west  dimensions  of  50-200  km,  often  extended  to  more  than 
1200  km  in  the  north-south  direction,  and  indeed  were  viewed  as  continuous 
irregularity  features  spanning  the  magnetic  equator  to  both  intertropical 
arcs.  The  depletions  appeared  after  sunset  and  were  observed  to  drift  to 
the  east  with  speeds  of  50-150  m/sec. 

O 

The  6300  A  airglow  imaging  technique  pioneered  by  Weber  et  al ,  (1978,1980), 
and  recently  extended  by  Moore  and  Weber  (1981)  to  include  images  at  7774  A, 
has  shown  that  the  equatorial  airglow  depletions  are  linked  unmistakably  to 


the  equatorial  F-region  bubbles  and  plumes  of  plasma  irregularities  that 
have  been  the  subject  of  intense  study  during  the  past  several  years 
(Matsushita  et  al . ,  1981).  Radio  and  rocket-borne  diagnostics  of 
equatorial  plasma  depletions  offer  very  detailed  "line  of  sight"  aspects 
of  the  phenomenon;  all-sky  optical  methods  offer  broad,  instantaneous 
coverage  of  the  features,  but  only  from  the  limited  altitude  range  where 
natural  airglow  occurs.  All  threemethods  (radio,  in-situ  and  optical)  are 
clearly  complementary.  The  imaging  technique  is  still  in  its  early 
phase,  and  promises  to  be  a  valuable  method  for  quantitative,  upper  atmos¬ 
pheric  science. 

In  this  paper  we  wish  to  describe  a  new  series  of  optical  measure¬ 
ments  obtained  from  ground-based,  photographic  methods.  The  results  suggest 
that  a  substantial  data  base  of  equatorial  depletion  characteristics  can  be 
obtained  during  a  relatively  brief  experimental  campaign,  and  that  the  air- 
glow  images  can  be  processed  to  yield  quantitative  descriptions  of  a  broad 
range  of  parameters  related  to  the  theory  and  practical  consequences  of 
equatorial  plasma  irregularities. 

In  the  following  Section,  our  all-sky  airglow  imaging  system  is  described 
In  Section  3,  we  describe  a  field  test  of  the  apparatus  during  an  equatorial 
observing  campaign  in  January  -  February  1981.  Observations  and  results  are 
discussed  in  Section  4,  and  a  summary  of  findings  appears  in  Section  5. 


2.  A  LOW-LIGHT-LEVEL  AIRGLOW 


IMAGING  SYSTEM 

In  equatorial  regions  the  normal  background  of  6300  A  radiation  is 
relatively  weak  (50-100  R).  In  order  to  detect  regions  of  depleted  airglow 
in  relatively  short  integration  times  (-'1  minute),  an  intensified  camera 
system  is  required.  The  optical  system  developed  for  this  purpose  is  shown 
in  schematic  form  in  Figure  1.  While  the  system  has  interchangeable  all-sky 

i 

(180°  field  of  view)  and  narrow  (60°  field  of  view)  lenses,  all  of  the  observa¬ 
tions  reported  on  here  were  taken  using  the  1 80 ^  field  of  view  option.  In  this 
mode,  a  standard  NIKKOR  8  mm  f/2.8  "fisheye"  lens  produces  a  circular  image  of 
the  sky  23  mm  in  diameter.  The  light  from  this  image  is  collimated  by  a  100  mm 

o 

diameter  f/2.8  lens  before  passing  through  a  6300  A  interference  filter  to  reject 

unwanted  background  radiation  from  the  Moon  and  stars.  The  system  was  designed 

0 

for  a  filter  having  a  very  narrow  bandwidth  (6  A  FWHP);  in  the  observations 

o 

reported  here,  an  alternate  12  A-wide  filter  was  used. 

After  passing  through  the  filter,  the  monochromatic  light  is  re-imaged  on 
the  photocathode  of  a  second  generation  image  intensifier  whose  luminous  gain 
is  approximately  60,000.  The  output  phosphor  is  then  imaged  on  35  mm  film 
using  two  NIKKOR  lenses  front-to-front  and  a  NIK0N-F2  camera  body.  All  of  the 
images  discussed  in  this  paper  were  taken  at  10  minute  intervals  using  a 
30-second  time  exposure. 


3.  THE  ASCENSION  ISLAND  EQUATORIAL 
IONOSPHERE  OBSERVING  CAMPAIGN  OF  1981 


The  Air  Force  Geophysics  Laboratory  organized  a  mul ti -techi nque  observing 
program  on  Ascension  Island  (7.95'  S,  14.38'W)  during  the  period  24  January 
to  10  February  1981.  Boston  I’niversity  ’  s  all-skv  airglow  imaging 
system  operated  from  27  January  to  8  February  and  obtained  useable  airglow 
data  from  11  of  the  13  nights.  Figure  2  shows  the  location  of  Ascension 
Island  with  respect  to  geographic  and  geomagnetic  coordinates.  It  can  be 
seen  that  Ascension  Island  sits  close  to  the  southern  crest  of  the  eguatorial 
anomaly,  essentially  beneath  the  southern  intertropical  arc.  The  two  circles 
indicate  the  coverage  of  the  camera  system  for  zenith  distances  of  75'  and 
90  ,  respectively.  In  Figure  3,  a  samole  photograph  is  shown  together  with 
three  types  of  overlay-grids  used  to  examine  the  image  features.  Note  that 
the  field  of  view  for  this  and  all  subsequent  photographs  exhibits  the  usual 
inversion  of  east-west  coordinates.  The  field  of  view  also  contains  a  group 
of  four  telephone-pole-size  antenna  poles,  with  connecting  wires.  These  poles 
are  reproduced  on  the  grids,  and  serve  as  convenient  orientation  markers,  as  do 
the  two  building  lights  on  the  eastern  he  izon.  The  widths  of  the  wires 
illustrate  the  resolution  achieved  by  the  system  (  1Y). 

The  grids  presented  in  Figure  3  contain  important  information  on  the 
spatial  extent  and  characteristics  of  the  images.  Note  that  the  zenith  angle 


grid  (lower  right)  describes  a  nearly  linear  relationship  between  zenith 
angle  and  image  size.  The  geographical  grid  above  it  shows  that  a  zenith 
angle  of  75  encompasses  approximately  -  8  of  latitude/longitude  from  the 
zenith,  which  is  equivalent  to  a  horizontal  diameter  of  -1800  km  thrc-uuh 
the  zenith  at  300  kir.  Between  zenith  angles  of  75  to  SO  ,  an  additional 

'8  of  latitude  longitude  are  encompassed,  but  coiiipression  effects  render  this 
region  difficult  to  work  with  in  any  quantitative  way.  It  is  important  tc 
note,  however,  that  the  geomagnetic  equator  falls  within  the  75'  -  90  region 
o4“  the  Meld  of  view,  as  illustrated  in  the  lower  left  grid. 

As  a  further  aid  to  interpreting  the  photographs ,  Figure  <3  gives  a 
schematic  representation  of  the  airglow  viewing  geometry  computed  for  the 
magnetic  meridian  of  Ascension  Island.  The  broad  band  near  a  height,  of  300  km: 
depicts  the  approximate  reqion  of  6300  A  emission  recorded  by  the  optical 
system.  From  north- to-south  horizon,  this  region  spans  3600km,  and  shows 
that  the  airglow  features  recorded  describe  characteristics  at  the  base  of 
many  field  lines.  Thus,  at  a  zenith  angle  of  85  to  the  south,  the  images 
recorded  from  airglow  characteri sties  near  300  km  pertain  to  field  lines 
that  cross  the  magnetic  equator  near  1850  km.  Images  recorded  near  the 
zenith  pertain  to  flux  tubes  crossing  the  equator  near  700  km,  while  images 
recorded  at  zenith  angles  -85°  to  the  north  arise  from  300  km  features  at 
the  magnetic  equator  itself.  Images  taken  from  such  a  geometry  have  the 
interesting  quality  that  features  seen  from  top-to-bottom  on  the  photographs 
have  a  one-to-one  correspond emu-  to  altitude,  bottom- to-top  (-300  km  to  1800  km), 
above  the  equator.  When  interpreted  in  this  way,  the  optical  imaging  technique 
can  be  seen  to  yield  information  about  the  height-dependence  of  plasma  depletions 
above  the  equator,  as  well  as  the  broad  morphology  aspects  inherent  in  an 
all-sky  photograph. 
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IV  R:  SUITS 

Figure  5  contains  a  summary  of  a  niqht's  observations  from  •  ■ rent  ion 
Island.  Located  close  to  15  local  time  at  Ascension  is  approximately  1 
nour  earlier  than  the  U.T.  values  shewn  in  the  figure.  The  bright  image 
of  the  setting  noon  is  seen  in  the  upper  right  corner  of  the  frames  at  21:0, 
and  21:10  U.T.,  and  the  bright  dot  near  the  two  building  lights  in  trie  lower 
left  results  from  the  combined  images  of  Jupiter  and  Saturn.  The  availability 
of  useable  airqlow  data  for  depletion  related  studies,  as  shown  here  for 
7-8  February  1981,  is  typical  for  that  obtained  on  other  nignts.  The  back¬ 
ground  airqlow  usually  increased  after  21 :00 U.T.  .  generally  revealing  well- 
formed  depletions.  After  01:00  U.T.,  the  depletions  gradually  disappeared 
and  observations  were  ended  for  the  night. 

Figure  5  exhibits  a  rich  display  of  equatorial  airqlow  depletion  shapes 
and  forms.  Frame  22:20  U.T.  alone  shows  narrow  bands,  twisted  features,  air- 
glow  islands  and  a  "wishbone'  formation.  Later  in  the  evening,  single  large 
depletions  dominate  the  scene.  Figure  6  gives  further  examples  of  narrow  depletion 
bands  (1  Feb '81,  top  panels)  and  a  second  "wishbone"  feature  (2  Feb '81,  bottom 
panels).  Figures  5  and  6  illustrate  several  morphology  features  that  can  be 
characteri zed  using  simple  overlay/grid  data  reduction  techniques.  These 
include  depletion  occurrence  rates,  alignment  characteristics  and  large-scale 
drift  speeds.  More  detailed  image  processing  methods  applied  to  special  features 


(such  as  the  ’wishbone")  and  to  individual  depletion  edne  effects  will  be 
described  in  a  following  section. 

4.1.  Depletion  Occurrence  Rates  and  Spatial  Coverage 

Durinq  the  13-night  optical  observing  program  on  Ascension  Island,  two 
nights  were  completely  cloudy  and  no  6300  A  emission  was  observed.  On  three  of 
the  remaining  eleven  nights,  there  was  a  bright  6300  A  airglow,  but  no  depleted 
regions  were  observed.  Thus,  on  eight  of  the  eleven  clear  nights  ('70  ),  well- 
formed  airglow  depletions  were  photographed  throughout  the  evening  hours. 

As  noted  in  the  previous  section,  regions  of  depleted  airglow  usually 
appeared  after  20:00  L.T.,  disappear  after  midnight,  and  occurred  most  abundantly 
during  the  20:30-23:30  L.T.  hours.  Baumgardner  and  Klobuchar  (1981)  analyzed 
each  photograph  containing  depletions  with  an  equal -sky-area  grid  in  order  to 
estimate  the  portion  of  the  sky  covered  by  depletions.  To  avoid  compression 
effects  near  the  edge  of  each  photograph  they  limited  their  analysis  to  zenith 
anqles  down  to  70°.  Their  results,  summarized  in  Figure  7,  show  that  on  several 
nights  the  depleted  regions  covered  over  30%  of  the  sky  for  extended  periods,  and 
on  a  few  occasions  over  40  of  the  sky.  Averaging  over  the  samples  available, 
one  concludes  that  on  nights  when  depletions  occur,  approximately  one-fourth  of 
the  sky  is  covered  by  depletions  for  at  least  two  hours. 

4.2.  Depletion  Alignments 

In  the  first  study  of  wide-angle  imaging  of  equatorial  airglow  depletions, 
Weber  et  al.,  (1978)  described  the  depletions  as  "north-south"  aligned.  Their 
observations  were  made  along  a  magnetic  meridian  near  Lima,  Peru,  from  an 
aircraft  that  kept  wi thin  ±3°  of  the  magnetic  equator.  They  noted  that  the 
depletions  were  basically  aligned  with  the  magnetic  meridians  close  to  the 
equator,  but  that  a  skewness  was  evident  for  depletion  ends  extending  far  from  the 
equator.  In  their  subsequent  study  (Weber  et  al . ,  1980),  the  aircraft  flew  in 
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a  region  above  and  to  the  north  of  Ascension  Island,  and  recorded  several  cases 
of  the  poleward  ends  of  depletions  curving  noticeably  to  the  west.  As  described 
in  Figure  4,  ground-based  photographs  from  Ascension  Island  are  particularly 
useful  for  documenting  alignments  since  the  magnetic  equator  appears  near  the  very 
top  of  the  photographs.  Hori zon-to-hori zon  depletions  in  the  N-S  direction 
thus  extend  over  3000  km  from  the  equator,  and  departures  from  magnetic  meridians 
can  be  seen  to  grow  from  top-to-bottom.  For  example,  consider  the  depletion 
shown  in  Figure  3.  The  central  depletion  begins  due  north  on  the  magnetic 
equator.  The  grid  below  it  shows  that  if  the  depletion  were  aligned  along  a 
magnetic  meridian,  it  would  curve  to  the  left  in  the  photograph.  Instead,  as 
one  progresses  down  the  axis  of  the  depletion,  it  appears  on  meridians  2°,  then 
4°,  and  at  the  bottom,  up  to  8°-10u  to  the  west.  Since  features  at  the  bottom 
of  the  photograph  relate  to  altitudes  far  above  the  equator,  this  westward  tilt 
must  be  the  optical  manifestation  of  the  westward  tilts  of  plasma  depletion 
plumes  recorded  by  incoherent  scatter  radar  (Woodman  and  La  Hoz,  1976;  Tsunoda, 
1980)  and  by  in-situ  probes  (McClure  et  al,.  1977).  The  end  result  of  these 
strong  westward  tilts  is  that  from  a  site  such  as  Ascension  Island,  where  the 
magnetic  declination  is  large  ('20°)  and  to  the  west,  the  depletions  appear  to 
be  geographical ly-al i gned  in  the  North-South  direction.  This  is,  of  course,  an 
artifact  of  the  geometry,  and  not  a  physically  important  characteristic. 

4.3.  Large-Scale  Plasma  Drifts 

Weber  et  al . ,  (1978)  described  several  cases  of  airglow  depletions  that  were 
observed  to  drift  eastward,  with  speeds  of  50-100  m/sec.,  for  up  to  several  hours. 
Figures  5  and  6  show  that  eastward  drifts  are  clearly  visible  in  10-minute  spaced 
photographs  from  Ascension  Island.  Using  the  eight  nights  of  airglow  images  that 
contained  depletions,  we  computed  eastward  drift  speeds  for  optical  features 
crossing  the  zenith  as  a  function  of  local  time.  Using  photographs  taken  10 


minutes  before,  on  ,  and  after  the  hour,  drifts  were  computed  for  features 
(dark  or  bright)  close  to  the  zenith,  usinq  an  assumed  height  of  300  km.  In 
some  cases,  when  a  photograph  was  rich  in  detail,  two  sets  of  speeds  were 
computed  for  a  given  hour.  The  results  of  this  analysis  are  given  in  Figure 
8,  where  the  average  speeds  (with  numbers  of  observations  in  parentheses)  are 
plotted  as  a  function  of  local  time.  Figure  8  also  contains  several  sources 
of  comparison  with  measurements  made  by  other  techniques: 

(1)  The  Sipler  and  Biondi  (1978)  results  refer  to  the  average  thermospheric 
wind  detected  by  the  Fabry-Perot  i nterferometer  method,  from  five  days  of 
observations  at  Kwajalein,  durinq  Auqust-September  1977.  The  Sipler  et  al . 

(1980)  results  refer  to  similar  measurements  during  five  days  in  July  1979. 

(2)  The  Aarons  and  Whitney  (1980)  results  were  obtained  by  analyzing  spaced 
receiver  scintillation  records  from  13  nights  in  March  1977  at  Ancon,  Peru. 

(3)  The  Fejer  et  al .  (1981)  results  were  obtained  from  incoherent  scatter 
observations  from  300-400  km,  during  several  years  of  operation  at  Jicamarca. 
Peru. 

The  results  of  Figure  8  come  from  a  variety  of  methods,  locations  and 
solar  flux  conditions,  and  contain  a  fair  degree  of  experimental  uncertainties. 
Nevertheless,  there  is  an  encouraging  agreement  between  the  local  time  dependence 
and  magnitudes  of  the  drifts  obta.ned  from  the  airglow  imaging  method  and  the 
other  traditional  approaches  to  F-region  dynamics.  Theoretical  considerations 
suggest  that  nighttime  plasma  drifts  near  the  equator  should  become  approximately 
equal  to  the  F-region  neutral  winds  which  generate  them  via  electrodynamic 
dynamo  action  (Rishbeth,  1971,  1981;  Heelis  et  al.,  1974).  The  pattern  of  east¬ 
ward  speeds  summarized  in  Figure  8  supports  that  view. 

A  close  examination  of  well-defined  features  in  Figures  5  and  6  (e.g., 


the  wishbone  formations)  reveals  that  northward  drift  speeds  are  also  evident 
from  the  airglow  images.  These  drifts  are  somewhat  more  difficult  to  extract 
using  overlay/grid  techniques  due  to  the  inherent  difficulty  in  trying  to 
identify  small  displacements  in  nebulous  features.  A  few  clearly-identified 
cases  were  obtained,  however,  and  they  yielded  northward  speeds  of  =50  m/sec 
from  22:00  to  23:00  L.T.  During  this  time  period,  eastward  speeds  average 
=130  m/sec,  and  thus  a  total  horizontal  ExB  drift  of  =140  m/sec  is  consistent 
with  the  declination  of  approximately  20"  near  Ascension  Island. 

4.4  Image  Processing  Techniques 

The  quality  of  photographs  capable  of  beinq  obtained  from  all-sky  airglow 
imaging  sugqests  that  computer  processed  photometric  methods  can  be  used  to 
extract  detailed  quantitative  information  from  the  photographs.  In  order  to 
test  these  procedures  with  our  initial  field  data,  the  density  distributions  on 
our  35  mm  imaqes  were  accessed  using  a  two-dimensional  scanning  E I K0N I X  micro¬ 
densitometer.  The  resulting  digital  output,  formated  to  a  128x128  pixel  array, 
were  converted  to  an  intensity  plane  using  well-known  film  calibration  procedures. 
Full  details  of  the  method  have  been  described  in  Chapter  2. 

In  order  to  demonstrate  the  types  of  digital  processing  capable  with  a 
relatively  simple  system,  each  10-minute  frame  from  two  nights  of  observations 
(1-2  and  2-3  February  1981)  were  digitized.  The  numerical  average  of  all  those 
frames  provided  a  simple  way  to  characterize  the  average  influences  of  the  Van 
Rhijn  and  vignetting  effects,  plus  field  contaminations,  such  as  building  lights 
and  telephone  poles  and  wires.  The  subtraction  of  a  single  frame  from  the 

average  then  yields  a  "difference  oicture"  which,  when  added  to  a  suitable 
numerical  offset,  provides  a  very  convenient  wa.'-  to  characterize  the  fundamental 
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aspects  of  airglow  depletions.  Digital  arrays  such  as  these  can  then  be 
used  to  construct  three-dimensional  (x,y,  intensity)  airglow  relief  maps, and 
associated  two-dimensional  contour  plots.  Figures  9  and  10  offer  examples  of 
the  technique.  One  can  see  that  the  method  captures  the  essence  of  a  photo¬ 
graph,  and  provides  the  added  benefit  of  quantitative  intensity  variations 
free  from  geometrical,  lens  and  field-of-view  effects.  The  contours  in 
Figure  9,  for  example,  show  that  the  depths  of  the  depletions  decrease  from 
west  to  east,  while  in  Figure  10,  the  wishbone  depletion  is  seen  as  a  continuous, 
bifurcating  formation  extending  several  thousand  kilometers  from  the  geo¬ 
magnetic  equator.  The  gr'.ds  in  Figure  3  show  that  the  left  or  eastern  branch 
of  the  wishbone  is  nearly  parallel  to  a  magnetic  meridian,  while  the  western 
branch  appears  to  be  in  the  geographic  north-south  plane,  as  discussed  earlier. 
Figure  11  contains  a  schematic  view  of  the  3-dimensional  characteristics  of 
field-aligned  depletions  that  could  account  for  such  an  "airglow  wishbone." 

The  digital  arrays  obtained  from  airglow  images,  such  as  in  Figure  9 
(1  February  1981  at  22:00  U.T.),  can  be  used  to  examine  the  temporal  evolution 
of  equatorial  depletions  in  a  high  spatial  resolution  format.  When  an  indivi¬ 
dual  frame  is  divided  by  the  normal i zed  composi  te  average  frame,  the  field  of 
view  is  corrected  for  van  Rhijn  and  vignetting  effects,  and  the  depletions 
appear  as  strong  features  on  a  relatively  uniform  background.  An  east-west 
(or  so-called  prime-verti cal )  scan  through  the  resultant  array  can  then  be 
used  to  characterize  the  instantaneous  cross  sections  of  all  the  features  in  the 
photograph.  Figure  12  contains  a  series  of  22  prime  vertical  scans,  spanning 
3’2  hours  of  observations  in  10-minute  steps,  from  the  night  of  1-2  February  1981. 
The  pixel  geometry  along  the  east-west  direction  has  been  transformed  into 
a  linear  format  in  kilometers  (at  h  =  300  km)  in  order  to  facilitate  comparisons 


of  features  photoqraphed  at  zenith  angles  from  0  to  75  east  and  west. 

Such  a  transformation  removes  the  image  compression  effects  that  occur 
towards  the  edge  of  a  photograph.  However,  since  the  pixel  spacing  was 
linear  in  the  digitization  process,  the  transformation  to  spatial  coordinates 
causes  resolution  to  deteriorate  badly  at  large  zenith  angles.  Figurel2 
therefore  exhibits  greatest  detail  within  ±600  km  of  the  zenith,  with  adequate 
resolution  for  broad  features  out  to  ±1200  km  (or  a  zenith  angle  of  75°). 

The  east-west  airglow  image  scans  in  Figure  12  are  plotted  with  successive 
75-unit  offsets  on  each  side  of  the  figure.  The  prime  vertical  scan  from  the 
first  photograph  of  the  evening  reveals  a  low  level  of  airglow  (-25  units)  that 
is  approximately  uniform  from  east  to  west.  As  the  F-region  descends,  back¬ 
ground  airglow  and  well-formed  depletions  emerge.  The  gradual  east-west  slope 
in  the  background  level  reflects  different  local  time  regions  across  a 
photograph  (i.e.,  it  is  latv.r  and  therefore  brighter  in  the  east).  Peak  back¬ 
ground  levels  occur  between  22:00-23:00  LT,  showing  approximately  a  7-fold 
increase  over  the  levels  recorded  earlier  in  the  evening.  As  the  F-region 
recombination  mechanisms  proceed  towards  midnight,  background  levels  decrease 
as  bottomside  plasma  is  consumed;  the  remnant  background  levels  now  remain 
high  towards  the  western  horizon  where  the  local  time  is  still  pre-midnight. 

As  seen  in  Figure  5,  the  prominant  airglow  depletions  tend  to  disappear  near 
midnight. 

Digitally  processed  image  scans,  as  shown  in  Figure  12,  can  be  used  to 
track  individual  depletions  over  their  3  to  3'2  hour  residence  time  within  the 
all-sky  field  of  view.  For  example,  the  depletion  that  appears  near  700  km  west 
of  Ascension  Island  at  20:40  LT  drifts  eastward  with  a  speed  of  :140m/sec  at 
21:00,  crosses  the  zenith  at  :120m/sec  an  hour  later,  and  finally  reaches  700  km 


cast  of  Ascension  Island  at  23:50  LT,  for  an  average  speed  of  '120  m/sec  from 
west  to  east. 

The  prime  vertical  scans  may  also  be  used  to  describe  the  cross-sectional 
characteristics  of  individual  ai rglow/pl asma  depletions.  For  example,  the 
scan  at  21:00  LT  in  Figure  12  depicts  three  depletions  that  are  2-2l2  times  less 
intense  than  the  background  level.  The  scan  at  23:00  LT  exhibits  similar 
factor  of  2  depletions  upon  a  background  that  is  now  twice  as  bright  as  the 
21:00  LT  scan.  Weber  et  al . ,  (1980)  reported  depletions  of  similar  magnitude 
using  an  airborne  photometer  to  monitor  zenith  intensities  of  the  drifting  deple¬ 
tions.  From  Figure  12,  it  becomes  clear  that  digital  processing  of  all-sky 
images  offers  an  attractive  method  of  separating  spatial  and  temporal  changes, 
while  still  preserving  the  quantitative  benefits  of  narrow  beam  photometer 
observations . 

The  east-west  walls  or  longitudinal  gradients  on  both  sides  of  the  deple¬ 
tions  are  particularly  well  defined  in  a  series  of  prime  vertical  scans.  In 

Figure  12,  over  50  depletion  cross-sections  are  shown  from  22  photographs  taken 

o 

over  a  3\  hour  period.  In  the  majority  of  cases,  the  6300  A  intensity  gradient 

on  the  western  wall  is  sharper  than  the  corresponding  gradient  on  the  eastern 

0 

wall.  Weber  et  al.,  (1980)  have  described  a  way  of  converting  6300  A  intensity 

to  the  average  electron  densi ty  within  the  airglow  emission  region,  and  thus  the 

brightness  gradients  in  Figure  12  should  be  good  indicators  of  the  actual  bottom- 

side  N  gradients  that  define  "plasma  bubbles."  The  steepness  history  of  a 
e 

western  wall  is  quite  evident  in  the  single,  long-lived  depletion  described 
above;  it  appears  on  the  western  horizon  in  Figure  12  shortly  before  21:00  LT 
transits  the  zenith  near  22:00  LT,  and  passes  to  the  eastern  horizon  after  23:00 
LT.  Throughout  this  period,  the  western  wall  is  steeper  than  the  eastern  wall. 
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It  should  be  noted,  however,  that  the  pronounced  steepening  of  both  walls  as 
the  depletion  passes  overhead  in  simply  an  effect  a  viewinq  geometry:  at 
large  zenith  angles,  the  slant,  integral  measure  of  any  trough-like  feature 
will  always  characterize  vertical  walls  as  gradual  boundaries  in  comparison  to 
overhead  observation  of  the  same  feature  (a  well-known  problem  encountered, 
for  example,  in  attempts  to  measure  the  edges  of  the  so-called  mid-latitude 
ionospheric  trough). 

A  final  point  to  be  mentioned  with  respect  to  Figure  12  concerns  the 
intensity  fluctuations  that  appear  in  the  bright  background  levels  near  the 
zenith  shortly  before  and  after  23:00  LT.  If  characterized  as  enhancements  upon 
a  background,  they  appear  as  20  increases  in  brightness,  a  value  that  is 
several  times  larqer  than  the  single  pixel -to-pi xel  noise  (-5  )  seen  in  tne 
briqht  reqi on  near  21 : 00  LT.  A  20  brightness  variation  is  still  too  small  to 
see  on  photographs  printed  to  portray  the  factors  of  2  to  3  difference  found 
in  depletion-to-background  effects.  The  enhancements  appear  to  be  tens  of 
kilometers  in  east-west  extent,  with  lifetimes  of  tens  of  minutes.  Further 
analysis  is  clearl  v  required  to  determine  if  they  are  actual  "fingers"  or  stri.it- ions 
of  enhanced  plasma  extending  from  the  sharp  western  wall  of  the  neighboring 
depletion,  evidence  of  Tin's  in  the  airglow  generating  region,  or  some  unknown 
problem  associated  with  experimental  or  digitizing  methods.  For  the  present, 
where  large-scale  depletion  cha rac ten s ti cs  are  our  main  concern,  it  is  best  to 
consider  these  fluctuations  as  noise  in  the  background  levels  that  outline  the 


depletions. 


5.  SUMMARY  AND  CONCLUSIONS 


A  new  low-light-level,  all -sky,  airglow  imaging  system  was  developed  to 
carry  out  the  first  ground-based  studies  of  the  optical  signatures  of  equatorial 
plasma  depletions.  The  basic  contribution  of  the  techniques  rests  on  the 
assumption  that  6300  A  emission  comes  from  a  relatively  narrow  altitude  range 
in  the  bottomside  F-reqion,  and  thus  a  region  of  depleted  airglow  represents  the 
base  of  a  flux  tube  with  severely  depleted  plasma  content.  A  narrow  band  of 
depleted  airglow  extending  away  from  the  equator  therefore  relates  to  a 
series  of  depleted  flux  tubes  that  cross  the  equator  at  increasingly  qreater 
heights. 

A  two-week  field  test  of  the  optical  system  designed  to  conduct  equatorial 
airqlow  depletion  studies  was  carried  ou^  on  Ascension  Island  in  early  1981. 

The  initial  results  to  come  from  these  observations  may  be  summarized  as  follows 

1.  Airqlow  depletions  were  observed  on  70  of  the  clear  nights.  Thev  appeared 
after  20:00  LT,  occurred  most  abundantly  during  the  20:30  -  23:30  LT 
period,  and  disappeared  after  midnight.  On  the  average,  they  covered 
approximately  one-fourth  of  the  sky  for  at  least  two  hours. 

2.  The  depletions  often  extended  from  the  maqnetic  equator  to  beyond  the 
crest  of  the  southern  i ntertropical  arc  (  3000  km).  Individual  depletions 
had  east-west  cross  sections  of  100  to  several  100  km. 

3.  The  airglow  depletions  were  aligned  in  a  variety  of  directions,  all  with 
approximate  north-south  ori enta ti ons .  In  most  cases,  depletions  were 
not  aligned  with  a  magnetic  meridian,  but  rather  showed  a  progressive 
skewness  to  the  west  at  further  distances  from  the  equator. 

4.  Several  cases  of  apparently  twisting,  overlapping  and  bifurcating  deple¬ 
tions  were  observed. 
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5.  The  depletions  drifted  to  the  east  with  speeds  that  decreased  from 
190  m/sec  at  21:00  LT  to  '80  m/sec  at  01:00  LT. 

b.  Airglow  intensities  within  depletions  were  typically  40-50  of  the 
background  intensities  that  defined  them. 

7.  The  western  wall  of  a  depletion  usually  had  a  steeper  gradient  than 
the  corresponding  eastern  wall. 

8.  Possible  evidence  exists  for  small  regions  ('10  to  50  km)  of  enhanced 
airglow  embedded  in  the  bright  background  areas  between  depletions. 

The  results  outlined  above  support  and  extend  the  conclusions  reached  in  the 
initial  series  of  equatorial  airglow  depletion  morphologies  obtained  by  weber 
et  al .  ,  (1978,  1980)  and  Moore  and  Weber  (1981).  The  airqlow  imagine  technioue 
may  also  be  used  for  high  resolution,  quantitative  investigations  of  depletion 
characteristics  by  using  relatively  simple  image  processing  methods.  This 
approach  is  most  useful  for  extracting  specific  spatial  and  temporal  parameters 
related  to  the  mechanisms  responsible  for  the  equatorial  depletion  effects. 

Zalesak  et  al  ,  (1982)  have  recently  presented  a  comprehensive  review  of 
their  numerical  simulation  studies  of  the  nonlinear  evolution  of  the  gravita¬ 
tional  ly-dri ven ,  collisinnal  Rayleigh-Taylor  instability  considered  responsible 
for  equatorial  r^asma  bubbles.  Their  simulation  results  contain  many  features 
reported  in  earlier  experimental  work,  especially  incoherent  scatter  observations 
from  Jicamarca  (Woodman  and  laHoz,  1976)  and  Kwajalein  (Tsunoda,  1980).  The  radar 
results  are  often  depicted  in  maps  of  meter-size  irregularity  characteri sties  as 
a  function  of  space  and  time.  The  simulations  are  not  aimed  at  these  small- 
scale  perturbations ,  but  rather  at  the  large-scale  (100's  km)  characteristics 
of  the  envelopes  or  boundaries  of  the  small  scale  effects.  The  optical  results 
presented  here  have  a  spatial  extent  and  scale  length  that  are  more  closely 
related  to  the  simulation  results  carried  out  to  date.  Several  large-scale 
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perturbation  results  computed  by  Zalesak  et  al.,  such  as  depletion  widths,  tilts, 
bifurcations  and  steep  west  wall  gradients  below  hmax,  are  all  evident  in  the 
data  set  described  here.  Further  detailed  comparisons  of  simulation  results  and 
airglow  imaging  data  may  prove  useful  in  understanding  the  complex  processes 
that  control  the  evening,  equatorial  ionosphere. 
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FIGURE  CAPTIONS 


Schematic  view  of  low-light-level  camera  system  developed  for 
equatorial  airglow  depletion  studies. 

Location  of  Ascension  Island  with  respect  to  geographic  and 
geomagnetic  coordinate  systems.  The  two  concentric  rings  mark 
the  ’’all-sky"  camera  coverage  for  zenith  distances  of  76  and 
901'',  respectively.  The  heavy  dashed  lines  mark  dip  latitudes 
of  ±15  ,  giving  the  approximate  locations  of  the  crests  of  the 

equatorial  (Appleton)  anomaly  and  their  associated  intertropical 
airglow  bands.  Several  equatorial  observing  sites  are  shown  in 
relation  to  Ascension  Island. 

Sample  photograph  showing  equatorial  airglow  depletions  embedded 
in  a  weak  background  of  6300  A  emission,  togetner  with  three  scale- 
size  grids  giving  qeographic,  geomaanetic  and  zenith  angle  para¬ 
meters.  Note  the  inverted  east-west  directions,  the  markers  for 
telephone  poles  in  the  field  of  view,  ana  3  zenith  angle  calibra¬ 
tions  on  each  of  the  grids. 

Magnetic  meridional  geometry  across  the  equatorial  region  near 
Ascension  Island,  and  its  relationship  to  airglow  features  recorded 
in  all-sky  photographs,  (see  text.) 

A  series  of  6300  A  photographs  taken  at  10-minute  intervals 
during  the  night  of  7-8  February  1981.  (see  text  for  discussion.' 

Examples  of  a  1  i - sky  images  of  6300  A  airglcw  features  recorded 
at  20-minute  intervals.  The  top  panel  (1  February  1981)  shows  the 
motion  of  three  narrow  depletions  towards  the  northeast.,  while  the 
lower  panel  (2  February  1981)  snows  an  example  of  a  "wishbone1'  air- 
qlow  depletion.  The  bright  backum  .mi  levels  highlight  the  telephone 
poles  and  wires  that  serve  as  convenient  orientation  markers;  two 
buildings  with  rooftop  lights  obscure  a  small  portion  of  the  eastern 
horizon.  Resolution  limits  may  be  estimated  by  the  widths  of  the 
wires  and  the  comoined  images  of  Jupiter  and  Saturn  : the  bright 
elongated  feature  that  moves  towards  tne  west,  close  to  the  equate1', 
from  the  building  lights  on  the  eastern  horizon.) 

Summary  of  airglow  depletion  extent  with  respect  to  the  visible 
airglow  from  Ascension  Island  (from  Baumgardner  and  ±1 obuchar ,  ]9oli. 

Average  oastu.mi  drift  of  airglow  features  versus  local  time, 
derived  from  all-skv  images  recorded  from  January  to  February  1981. 

The  number  of  individual  drifts  used  to  form  the  average  speeds  appear 
in  parentheses.  Eastward  speeds  measured  by  other  techniques  ami  at 
other  locations  are  shown  for  comparison.  ; see  text.) 

Example  of  digital  processing  results  for  tne  all-skv,  b3l)0  A 
airglow  image  recorded  at  22:00  LIT  on  1  February  1981.  The  three- 
dimensional  relief  map  uses  a  "difference  from  the  mean’  technique 
to  correct  for  lens  and  field  of  view  effects;  the  two-dimensional 


Figure  9. 


contour  of  the  airglow  map  shows  depletion  effects  with  contours 
of  relative  intensity  (see  text). 

Figure  10.  Example  of  digital  image  processing  upon  the  "wishbone"  forma¬ 
tion  recorded  in  the  6300  A  airglow  image  recorded  at  00:30  UT  on 
3  February  1981.  (see  Figure  6  and  text.) 

Fiqure  11.  Schematic  representation  showing  how  a  wishbone-shaped  equatorial 
airglow  depletion  maps  along  geomagnetic  field  lines  to  outline 
plasma  depletion  characteristics  above  the  equator. 

Figure  12.  A  series  of  east-west  ("prime-vertical")  scans  from  the  digitul 
arrays  created  from  the  22  photographs  of  6300  A  airgnu  images 
recorded  on  the  night  of  1  February  1981.  The  data  are  corrected 
of  vignetting  andean  Rhiin  effects,  and  converted  to  east-west 
distances  assuminq  a  mean  heiqht  of  300  km  for  the  airqlow  emission. 
The  scans  at  20:00  LT  (bottom  left)  and  22:10  (bottom  riqhti  are 
plotted  in  relative  intensity  units,  with  10  subsequent  scans  plotted 
with  successive  73  unit  offsets. 
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CHAPTER  r> 


OPTICAL  SIGNATURE  OF  AN  IONOSPHERIC  HOLE 


INTRODUCTION 

Rocket  exhaust  induced  effects  upon  the  ionosphere  are  a  form  of  labora¬ 
tory-in-space  experimentation  that  can  be  used  to  address  a  broad  range  of 
chemical  and  dynamical  processes  of  interest  to  aeronomic  and  space  plasma 
research  (Mendillo,  1981).  While  radio  diagnostic  systems  have  been  the  main 
technique  to  monitor  such  effects  (Mendillo  et  al,  1975,  1980),  optical 
methods  can  provide  a  complimentary  way  of  observing  disturbance 
morphologies  (Kofsky  et  al,  1980;  Stone  and  Zinn,  1980).  A  recent  rocket  launch 
provided  a  new  "experiment  of  opportunity"  to  test  the  utility  of  low  light 


level,  image  intensified  photography  as  a  diagnostic  tool  for  ionospheric 
modification  experiments. 


OBSERVATION/- 

The  NOAA-C  weather  satellite  was  launcned  at  10:53  U.T.  from  Vandenberg  Air 
Force  Base  (34.6°N,  120. 6°W)  on  23  June  1981.  The  Atlas-F  launch  vehicle  followed 
a  rarely  used  ascent  profile  that  called  for  an  engine  burn  well  into  the  topside 
F-region  (see  Figure  1) .  A  temporary  observing  station  was  established  near  Edwards 
Air  Force  Base  (35.0°N,  117. 8°W)  to  make  radio  beacon  and  optical  measurements  of 
the  anticipated  effects.  The  site  was  chosen  in  order  to  have  the  ray  path  of  the 
136  MHz  beacon  from  ATS-1  pass  through  the  rocket's  exhaust  plume  at  350  km.  The 
VHF  radio  beacon  was  used  to  measure  the  total  electron  content  (TEC)  along  the 
ray  path  by  the  standard  Faraday  rotation  technique  (Titheridge,  1972).  The 
site  was  aLso  selected  on  the  basis  of  clear  sky  conditions  in  order 
to  provide  the  first  opportunity  to  make  optical  measurements  along  the  same  path 
as  a  TEC-detected  hole. 

The  optical  system  consisted  of  a  photometer  with  a  1°  field  of  view  and  a 
low  light  level  imaging  system  with  interchangeable  all-sky  (180°)  and  narrow- 
field  (60°)  lenses.  The  imaging  system  records  intensified  output  on  35  mm  film 
using  a  conventional  single  lens  reflex  camera.  All  observations  were  made  using 
6300  8  filters  that  have  a  very  narrow  half-power  width  (6  8),  a  condition  made 
necessary  by  the  presence  of  the  third  quarter  moon.  Control  data  were  taken  at 
6200  8  and  the  photometer  records  were  reduced  to  absolute  values  using  a  calibra¬ 
ted  source.  The  narrow  field  data  reported  on  here  were  obtained  using  8  second 
exposures . 

RESULTS 

Figure  2  contains  a  summary  of  the  electron  content  and  photometer  measurements 
made  during  the  event.  The  Faraday  rotation  data  show  a  very  rapid  depletion  oi  the 
F-region  followed  by  a  burst  of  airglow  f rom  the  oxygen  line  at  6300  8.  Figure  1 
presents  a  sample  of  the  imaging  results  using  the  b0°  1  told  oi  view  lens.  Note 
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the  posi.ci.on  ot  the  photometer  monitoring  point,  the  AfS-1  point,  and  the  times 

of  the  photographs  in  relation  to  the  photometric  curve  in  Figure  2.  The  images 

in  Figure  3  show  that  a  "smoke-ring"  or  "horsesh.M.— sn.»pcd"  wave  01  airgiuw  expand¬ 
ed  with  the  cloud  of  exhaust  molecules  responsible  for  the  F-region  depletions. 

This  effect  is  described  in  more  quantitative  fashion  in  Figure  4.  Digital  pro¬ 
cessing  of  the  images  was  carried  out  by  scanning  a  512  x  512  array  on  an  EIKON- 
1XSCAN  7800  system.  i'hese  data  were  corrected  subsequently  for  lens  effects  and 
smoothed  to  the  128  x  128  pixel  array  contoured  in  Figure  4.  The  horseshoe  pat¬ 
tern  is  quite  evident  on  a  log  (Intensity)  format,  and  the  lowest  level  contour 

(log  I  =  1 . 5  -  100  R)  gives  a  good  representation  of  the  outermost  visible  edge 
on  the  photograph.  The  full  evolution  of  this  contour  level  provides  an  overall 
summary  of  the  spatial  and  temporal  development  of  the  airglow  pattern  (see  Figure  5)  . 

DISCUSSION 

Rocket  exhaust  clouds  typically  include  molecules  of  H2 ,  H20  and  C02 •  These 
gases  react  with  the  dominant  F-region  ion  (0+)  to  form  molecular  ions  (0H+,  H20+, 

j  _L. 

H30  ,  02  )  at  rates  100  to  1000  times  faster  than  the  normal  0  reactions  with 
atmospheric  02  and  N2.  The  end  result  is  that  an  atomic  ion  plasma  is  converted 
to  a  molecular  ion  plasma  that  can  recombine  quickly  with  electrons.  Since  the 
rocket  plume  diffuses  rapidly  through  the  F-region,  the  enhanced  recombination 
process  occurs  over  a  broad  altitude  range  and  a  columnar  electron  content  deple¬ 
tion  is  observed.  The  dissociative  recombination  reactions  that  cause  the  hole 
can  leave  atomic  oxygen  atoms  in  excited  states  and  thus  to  a  burst  of  airglow. 
Theoretical  models  for  hole-making/airglow  scenarios  have  been  described  by 
Anderson  and  Bernhardt  (1978)  and  Zinn  and  Sutherland  (1980). 

Measurements  reported  here  of  the  evolving  airglow  cloud  can  be  used  to  test 
various  chemical,  diffusive  and  optical  parameters  related  to  plasma  modification 
mechanisms.  Two  aspects  are  addressed  in  this  preliminary  analysis:  diffusion 


and  airglow  production. 
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1.  Diffusion.  The  diffusion  of  modification  related  gases  in  the 

F-region  has  been  treated  by  Bernhardt  et  al.  (1975),  Mendillo  and 

Forbes  (1978),  Schunk  (1978),  and  Bernhardt  (1979).  Since  the  injected 

molecules  react  quickly  with  the  ambient  ions,  and  the  6300  R  emission  occurs 

with  the  subsequent  dissociative  recombinations,  the  expanding  airglow  pattern 

should  be  related  to  time  constants  associated  with  molecular  diffusion.  Using 

a  simple  diffusion  model,  one  can  show  that  the  time  of  peak  0+  loss  (t  )  is 

r  max 

related  to  distance  (  r  )  from  an  initial  source  as  t  =  rz/6U, 

max 

where  D  is  the  diffusion  coefficient  for  the  injected  molecules  (Mendillo  et  al . , 
1975)  . 

The  composite  contours  in  Figure  5  are  not  spherical  due  to  the  horizontal 
motion  of  the  rocket  and  the  exponential  nature  of  the  background  atmosphere. 
Nevertheless,  they  can  be  used  to  set  an  upper  limit  for  diffusion  using 
equivalent  spherical  arguments.  Figure  6  gives  the  effective  linear  expansion 
rates  along  the  major  and  minor  axes  of  the  contours,  together  with  a 
theoretical  equivalent  spherical  growth  rate  using  D  =  30  x  1.0  10  cm*/sec. 

This  value  is  consistent  with  the  diffusion  of  the  lightest  species  (U2)  through 
background  oxygen  at  350  km  and  T  =  1000°K. 

2.  Relationship  between  hole  formation  and  6300  R  airglow  emission. 

The  exhaust  gases  of  an  Atlas-F  rocket  include  the  following  molecular  injec¬ 
tion  rates  (in  units  of  1025  molecules/sec):  OH  =  1.5,  Ho  =  3.6,  COo  =  4.0  and 
1120  =  10.2.  Current  models  of  hole  formation  show  that  the  production 

of  the  0(1D)  state  of  oxygen  is  related  directly  to  the  dissociative 

+  +  + 

recombination  of  molecular  ions  (OH  ,  Oi  ,  HoO  )  formed  by  these  molecules.  The 
emission  of  6300  R  photons  is  described  by  the  Kinstein  coefficient  Aj  for  the 
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*D  ■+  3P2  transition  and  the  column  content  of  the  lD  state,  [0*] 


T. 


ET  -  Ai  j  [0*]  dh 


Ai [0*] T  photons 
cm2  sec 


At  x  10-9[0*]t  kR 


(1) 


For  the  unperturbed  ionosphere,  steady-state  conditions  are  usually  invoked  and 
[0*]T  is  computed  by  first  balancing  production  and  loss  at  each  height  (see, 
e.g.,  Forbes,  1970).  Since  a  modification  scenario  is  inherently  non-ste  dy- 
state,  such  an  assumption  may  not  be  used  here;  it  is  possible,  however,  to  use 
height- integrated  quantities  to  make  a  simple  estimate  of  the  time  dependence  of 
[ 0*1 T -  Thus,  if  motions  are  neglected, 

d[0*]T  -  P  -L  -  -k  dNr  -  (Aj,  +  QT)  [0*]r  (2) 

dt  dt 

The  first  term  represents  the  column  production  rate  of  0*,  where  k  is  the 
column-averaged  fraction  of  plasma  recombinations  that  yield  the  excited  0(*D) 
state, and  dNT/dt  is  the  observed  rate  of  change  of  the  columnar  electron  content. 
The  loss  term  derives  from  the  total  *D  emission  rate,  A^,  (*D  3Po,  3Pi,  3P2) , 
and  is  the  effective  total  columnar  quenching  rate.  Since  ANT  was  observed  to 
decrease  exponentially  with  an  effective  loss  coefficient,  0,  equation  (2)  may 
be  written  as 

,  k0ANT  e~St  -  (Aj.  +  Qt)  [0*Jt 
The  solution  to  (3),  when  coupled  to  (1),  yields 
_  k/ANj  X_10~1  re  -^.e'(AT  +  V'l 


8> 


The  coefficients  Aj  and  Aj,  are  fixed  at  0.0069  sec-1  and  0.0091  sec-1,  while 
ANt  and  3  are  obtained  from  the  observations  in  Figure  2:  g  =  0.0213  sec-1 
and  AN^,  -  16.8  10 12  el/cm2.  The  average  branching  ratio  k  for  the  number  of 
0*  per  recombination  and  the  column  average  quenching  coefficient  can  assume 
a  range  of  values  in  equation  (4) .  Three  pairs  of  (k-QT>  solutions  to  (4)  are 
plotted  in  Figure  7,  with  e.ach  case  chosen  to  yield  the  observed  peak  6300  % 
emission  corrected  for  atmospheric  extinction: 


E_  (t  )  =  8.9  kR  ate  =  ±  Ql>M. 

T  max  max  (Aj,  +  QT  -6) 


(5) 


It  becomes  clear  from  Figure  7  that  the  best  agreement  between  the  observed  and 

computed  8-(t)  curves  occurs  for  the  no-quenching  case.  With  ()„  =  0,  t  from 
r  T  max 

(5)  occurs  70  seconds  after  the  initiation  of  the  AN^  decrease;  the  observations 
in  Figure  2  show  that  the  peak  emission  occurred  nearly  120  seconds  after  the 
onset  of  the  hole.  Quenching  of  0*  by  molecular  species  is  of  the  form 
Q,j,  =  y[XY] ,  where  [XY]  is  the  average  molecular  concentration  within  the  hole 
and  y  is  approximately  5  x  10-11  cm3/sec  (Forbes,  1970;  Roble  et  al.,  1976). 

For  a  1000  °K  neutral  atmosphere  (CIRA,  1972),  N2  at  350  km  is  -  2  x  107  cm-3 
and  thus  -  0.001  sec-1or.only  10%  of  A^  in  equation  (5).  The  exhaust  molecules 
themselves  could  be  effective  quenchers,  especially  H.,0  (Hudson  and  Reed,  1979). 
This  could  be  important  at  early  times,  when  concentrations  v  107  cm- 3  are 
required  to  account  for  the  0  derived  above.  Nevertheless,  the  condition  Q.„  =  0 


in  equation  (4)  may  well  be  a  reasonable  first  approximation  for  estimating  the 
6300  A  airgiow  to  come  from  an  ionospheric  hole  formed  in  the  upper  F-region. 

Under  such  conditions,  approximately  14%  of  the  recombinations  yield  0 ( 1 D )  atoms. 
For  significantly  higher  quenching  rates,  proportionately  larger  \i,ldsof  0*  per 
recombination  are  required.  Such  adjustments  to  equation  (4)  cause  the  peak  emis¬ 
sion  to  occur  at  yet  earlier  times  and  tin*  decay  to  ambient  levels  occurs  much 


! aster  than  observed. 
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The  delayed  t  and  slow  return  to  ambient  6300  X  levels  shown  in  Figure  2 
max 

must  result  from  processes  not  considered  in  the  simple  formulation  leading  to 
equation  (4).  The  110  second  lifetime  for  0*  is  sufficiently  long  to  allow  for 
[0*]  to  change  by  diffusion  into  and  out  of  the  line  of  sight  monitored  by  the 
photometer.  Moreover,  the  diffusion  of  ambient  0+  into  the  molecular  rich  hole 
at  later  times  could  produce  additional  0*  that  are  not  detected  by  the  Faraday 
technique  (note  that  dN^/dt  -  0  from  11:00  -  11:15  U.T.  in  Figure  2,  and  thus 
further  plasma  recombinations  could  be  masked  by  a  matching  plasma  influx.)  In 
relating  the  production  of  [0*]  to  -k  (N^,/dt)  in  equation  (2),  there  is  an 
assumption  that  each  recombination  has  the  potential  of  producing  an  0*;  Anderson 
and  Bernhardt  (1978)  have  shown  that  in  a  dense  environment  of  H2  (as  occurs  at 
early  times),  molecular  ions  such  as  H20+  are  formed  that  recombine  very  rapidly 
(H30+  +  e~  ■+  H2O  +  H)  but  do  not  yield  0(JD)  states.  To  the  extent  that  this 
occurs,  the  production  of  (0*]^  in  equation  (2)  is  overestimated  at  early  times, 
and  thus  a  correction  in  the  proper  sense  is  obtained. 

The  assumption  of  a  high  efficiency  (k  =  1-1.3)  for  0(ZD)  production  per 
recombination  of  02+  (Roble  et  al.,  1976;  Hays  et  al.,  1978)  seems  at  odds  with 
the  low  column-average  k  values  depicted  in  Figure  7.  This  suggests  that  k  is 
a  tunction  or  relative  composition,  and  therefore  of  time  and  space.  Since  the 
yield  of  0*  is  high  from  02+,  the  CO?  component  of  the  plume  must  have  been  a 
minor  contributor  to  the  observed  TEC  decrease.  In  comparison  to  CO?,  the  other 
exnaust  gases  have  faster  reaction  rates  with  0+,  greater  diffusion  coefficients 
and  a  combined  abundance  'v  3-s  times  the  CO2  in  the  plume.  Thus,  while  H2  and 
H2O  are  very  efficient  hole-makers,  they  must  produce  0*  at  relatively  low 
efficiencies,  while  the  opposite  case  seems  to  occur  for  CO2. 

A  more  extensive  examination  of  this  and  a  similar  Atlas-F  event  using 
networks  of  TEC  and  6300  X  observations  will  be  described  in  a  1  it  i  r 
;.  i;i,  r.  The  single-site  data  presented  here  were  used  to  illustrate  the 
potential  benefits  to  be  gained  by  active  space  plasma  experiments  in  the 
iuiiu  plieric  domain.  I  ienii  -  ,  ■nt  .1  It:  ,  iiinmarv  :ew  >  •  -urli  work  . 
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